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ABSTRACT 



Gravitational lensing by massive galaxy clusters allows study of the population of intrinsically faint infrared galaxies that lie below 
the sensitivity and confusion limits of current infrared and submillimeter telescopes. We present ultra-deep PACS 100 and 160 /jm 
observations toward the cluster lens Abell 2218 to penetrate the Herschel confusion limit. We derive source counts down to a flux 
density of 1 mJy at 100 /im and 2mJy at 160 yum, aided by strong gravitational lensing. At these levels, source densities are 20 and 10 
beams/source in the two bands, approaching source density confusion at 160yum. The slope of the counts below the turnover of the 
Euclidean-normalized differential curve is constrained in both bands and is consistent with most of the recent backwards evolutionary 
models. By integrating number counts over the flux range accessed by Abell 2218 lensing (0.94-35mJy at 100 yum and 1.47-35mJy 
at 160;um), we retrieve a cosmic infrared background surface brightness of ~8.0and ~9.9 nWm"- sr"', in the respective bands. These 
values correspond to 55 ± 24% and 77 ± 31% of DIRBE direct measurements. Combining Abell 2218 results with wider/shallower 
fields, these figures increase to 62 ± 25% and 88 ± 32% CIB total fractions, resolved at 100 and 160 /jm, disregarding the high 
uncertainties of DIRBE absolute values. 

Key words. Surveys - Infrared: galaxies - Galaxies: evolution - Galaxies: high-redshift - clusters: general - Gravitational lensing: 
strong 
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1. Introduction 



The discovery of the cosmic infrared background (CIB) ( Puget et 
[HT 1996 Fixsen et al. 1996';'Lagache et al. 1999] ) has opened new 
perspectives on galaxy formation and evolution. A large number 
of the sources contributing to this far-infrared (FIR) CIB have 
been resolved in the mid-infrared (mid-IR) at 15//m with the 
ISOCAM instrument on ISO ( |Genzel & Cesarsky |2 000; Elbaz 



etal . 2002 ), and later at 24//m with MIPS on Spitze r ( |Fapovich| 



et al. 



,2004^ ^Dole et al. 



2504{ [Frayer et alT] [2006) 1. A striking 



result concerns the evolution of the infrared and submillime- 
ter (sub-mm) galaxy population: the infrared source counts are 
higher than no-evolution or moderate-evolution models and pro- 
vide strong constraints on the evolution of the bolometric energy 
output from distant galaxy populations. 
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* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



However, at FIR and sub-mm wavelengths, a much lower 
fraction has been resolved so far, because of the small aperture 
of telescopes, the prohibitive confusion limits, and the low sen- 
sitivity of available instruments. Nevertheless, stacking results 
(Dole et al. 2006) show that most or all of the FIR background 
is due to known high-z IR galaxies. 

The Herschel Space Observatory's (Pilbratt e t al.|[20T0) l 
compact point-spread-function (PSF) enables probes of the FIR 
emission of large samples of galaxies near their spectral en- 
ergy distribution (SED) peak, over a wide z range. The sensitive 
Photometer Array Camera &- Spectrometer (PACS) (Poglitschl 
et al. 2010) achieves imaging surveys of unrivalled depth at 70, 
100, and 160 yum in the PEP (PACS Extragalactic Probe) Key 
Program (PI D. Lutz). 

As pioneered by Sm ail et al.|( !1997) in the sub-mm, and later 
in the mid-IR (Altieri et al.||1999, Metcalfe et al.|[2003 ]), we 
are conducting a survey toward massive clusters as part of PEP, 
to resolve the FIR beyond the field-survey limit and penetrate 
the Herschel confusion limit with powerful gravitational lenses. 
This letter reports initial observations on the cluster Abell 2218. 
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A2218 Herschel PACS 160 micron 
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Fig. 1. (LEFT) PACS 100 ;um map toward Abell 2218 with lensed and foreground sources marked by blue circles. Other sources are 
identified cluster members. The rms contours at 0.7, 1.0, 2.0, and 4.0 mJy are overlaid in black. (RIGHT) PACS 160 jum map. The 
rms contours at 1.4, 2.0, and 4.0 mJy are overlaid in black. 



2. Observations 

The rich and massive lensing cluster Abell 2218 was tar- 
geted in the PEP science demonstration phase. This field has 
been intensively studied, including deep archival Hubble Space 
Telescope (HST)/ ACS and Spitzer (IRAC & MIPS) observa- 
tions, which allow simple and direct identification of the sources 
detected with Herschel/PACS. 

Our 13 hours of observations, centered on the cD galaxy 
at RA=16h35m51.84s, Dec=62dl4ml9.2s (J2000), used scan 
mapping (20"/s scan speed; 4' scan leg length; 20" cross-scan 
steps) over an area of about 6'x6', with coverage strongly peaked 
toward the center Though affected by relatively high overheads 
(reduced later in the mission), they are the deepest observations 
so far at 100 and 160 ;um. 

3. Data analysis 



The data were reduced using HIPE ( |Ott et al.|20T0| v2.0. build 
1328. The 10 Hz data cubes were processed with the standard 
PACS pipeline, along with custom procedures: 1) to remove in- 
terference patterns, transients from calibration blocks and track- 
ing anomalies, 2) to perform re-centering of positional offsets. 
Moreover, to remove detector drifts and 1// noise, a sliding 
high-pass filter was run on the pixel timelines with an iterative 
masking of the brighter sources. For more information, refer to 
Berta et al.^(2010) . Final maps are displayed in Figure[T| 

Source extraction was performed with the StarFinder PSF- 
fitting code (Diola iti et al.|2000| , down to the 3cr level. The total 
number of sources is reported in Table [T] Based on random ex- 
tractions from the residual images, the averaged Icr noise levels 
in the maps are 0.92 mJy at 100 /im and 1.61 mJy at 160 yum. 

We inspected all sources in the maps by eye. The source cata- 
log was cross-correlated with various redshift surveys for Abell 
2218 (LeBorgne et al.||19921 [Z Tegler et al.||200T] [Ebbels et al. 
1998 i Metcalfe et al. 2003 ; Kneib, priv. comm.) and with photo- 
metric redshifts of the Abell 2218 MIPS 24//m sources computed 



using HST/ ACS images, near-IR images and the 4 Spitzer/IRAC 
bands ( |Ciesla eFal. 200 9). All PACS sources were detected in 
the SpitzerfMIPS observation. This cross-correlation allowed us 
to classify sources as foreground, cluster, or background galax- 
ies. More than 90% of background sources with confirmed spec- 
troscopic or photometric redshift have Fif,o/Fioo Z 1, while 11 
out of the 13 confirmed cluster galaxies detected both at 100 and 
160 jum have Fi6o/^ioo < 1 . This allows classification of sources 
with Fi6o/^ioo ^ 1, but without redshift, as lensed-background 
galaxies. This is as also supported by their extremely faint mag- 
nitudes in the optical and/or their disturbed morphologies. 

Less than 30% (20%) of the sources at 100 //m (160 jum) are 
classified as cluster galaxies. This situation is very similar to the 
mi d-IR, wher e e.g. 'Alti eri eFaT] ([T999); [Metcalfe et al.| ( [2003) l; 
iEgami et al.| (2006) ; Hopwood et al.| ( |2010| l, showed that the 
large majority of far-IR sources are background, hence lensed, 
sources. The cluster core is virtually transparent in the FIR and 
acts as a natural telescope to provide a magnified view of the 
background sky. This increases the sensitivity and at the same 
time reduces the effects of source confusion; in particular, the 
cD galaxy is undetected in both bands. 

Four of the 7 sub-mm sources in the ultra-deep map of Abell 
2218 ( Knudsen et al.|2006) l are detected in the 160 /vm map, in- 
cluding the triple z=2.516 sub-mm source SMMJ16359H-6612 



Table 1. Properties of the A2218 PACS sources detection. 



mOfim leOyum 



PSF FHWM (arcsec) 



7.5 



11.0 



Icr average sensitivity (mJy) 
Icr central sensitivity (mJy) 
number of detected sources (> 3<t) 
number of lensed sources (> 3cr) 



0.92 
0.6 


1.61 

1.2 


98 
70 


94 

78 



lensed sources with spectroscopic redshift 9 



13 



lensed sources with photometric redshift 20 



17 
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( |Kneib et al.||2004a| l, the sources SMMJ163555.2-H 661150 at 
z= 1.034, and SMMJ163541.2-H661144atz = 3.1824 ("Knudsen 



et al. 2009 1. However, the highest redshift sub-mm source 



(z=4.048) SMMJ 163555. 5 +66 1300 and the two faintest sub- 
mm sources have not been detected, which is also true for the 
z=5.56 source of Ellis et al. ( 2001[ l and the z ~ 7 triple source of 
|Kneibetal.| ( r2004b| l. 

4. Lensing inversion and source counts 

Lensing acts in two ways on the background sky: i) it amplifies 
source brightness, typically by a factor of 2, but by as much as 
10 near critical lines; ii) it magnifies the area probed - and both 
the flux amplification and the space magnification are stronger 
toward the cluster core and increase with source-plane redshift. 



We exploited the detailed mass model of Abell 22 1 8 ( Kneib et 
3T]|1996 2004b I, which considers 8 multiply imaged systems, 
of which 7 have spectroscopically confirmed redshifts (among 
these the several high-z multiply-lensed sources mentioned at 
the end of Section 3). 

We used the publicly available software, LensTool, described 
in Julio et al. ( 2007| l to compute the source amplification factors 
and the surface magnification due to lensing. Where redshifts 
were not available for sources judged by morphology or intrin- 
sic faintness to be in the background, a value of z = 0.7 was 
assigned, corresponding to the median redshift of field sources 
having spectroscopic or photometric redshift. This has little im- 
pact on the lensing amplification factor and area magnification, 
as illustrated in Figure l2] This figure shows the area mapped to 
a given sensitivity in the image plane (i.e. on the apparent sky), 
compared to the corresponding areas in source planes at redshifts 
0.7, 1, 1.5, and 2.5. For redshifts above 0.7, the lensing surface 
magnification depends only weakly on redshift. The lensing am- 
plification factor for a given source on the other hand may be 
more sensitive to the redshift. However, setting all (morpholog- 
ically classified) background sources without known redshift to 
z=l (instead of z=0.7) did not change the shape of the counts 
significantly, because it shifted only a few sources among adja- 
cent flux bins. 

Figure [2] expresses the fact that some small, highly-lensed 
regions of the source plane are mapped onto larger areas in the 
image plane (the apparent sky), with the result that, over a few 
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Fig. 2. The area of the source plane at lOOyum, as a function of 
the 3cr sensitivity, for redshifts z=0.7, 1, 1.5, and 2.5, compared 
to the corresponding area in the image plane (i.e. apparent sky - 
the solid line). 



arcmin , flux densities in the range 1 to 3 mJy become accessible 
by virtue of the lensing effect. 

Ten sources (12 apparent sources, including all images of 
SMM J16359+6612) have lens-corrected fluxes below 3 mJy at 
lOOjum and 13 (15 apparent) below 5.7 mJy at 160;um, with 
these quoted limits the 3cr sensitivities achieved on the GOODS- 
N field (Bert a et al.||201()l l. Such sources would most likely 
not be detected even in upcoming deeper blank-field surveys, 
like GOODS-S (deeper by a factor of 2). The triple source 
SMM J16359-I-6612, for instance, is ampUfied by a factor of 45, 
in total. It was counted only once, as number counts refer to the 
source plane. 

By correcting for lensing amplification, surface magnifica- 
tion eff'ects, contamination by cluster galaxies, and non-uniform 
sensitivity of our maps, we can derive number counts at 100 //m 
and 160 /im. Because of the non-uniform sensitivity of the maps 
on the sky and the lensing eff'ect, different areas on the sky are 
surveyed to diff'erent depths in the source plane. The object den- 
sity per flux bin was computed using gain-dependent surface ar- 
eas. Incompleteness affects the counts progressively below ap- 
parent fluxes of 6mJy at lOOyum and 9mJy at 160yum. We re- 
stricted these counts to 4cr detections where the completeness of 
measurement is typically 80%. This avoids potentially compli- 
cated folding of completeness correction with lensing correction, 
not justified by the relatively small numbers of sources and rel- 
atively large statistical error bars associated with these sample 
counts on a single lensing cluster. A comprehensive complete- 
ness analysis will be required when combining the observations 
of a sample of massive cluster lenses to extract the galaxy num- 
ber counts below the 1 mJy level. 

Our counts are barely affected by confusion at 100/vm, with 
a source density of 20 beams/source in the image plane (for 
the Lagache et al.|200"3 definition of the beam). At 160;um the 
density is as high as 10 beams/source, with even fewer beams 
per source in the central area of high-lensing. Hence, the Abell 
2218 catalog is affected by source confusion at 160 /^m (16.7 



beams/source from Dole et al. ( 2003 i) and the high density of 
detected sources prevents the extraction of fainter objects. 

The source counts, corrected for cluster contamination and 
lensing effects, in both the 100/im and 160jum bands, are pre- 
sented in Figure[3| normalized to the Euclidean slope {dNIdS oc 
S^^'^). Error bars consider Poisson statistics only, as flux uncer- 
tainties are minimal by comparison. 

5. Discussion 

The Abell 2218 HerschelfPACS maps are the deepest FIR maps 
to date. Covering an area of 6 x 6 arcmin^ over the cluster 
core, we have surveyed the region where strongly lensed back- 
ground sources are present. At fluxes higher than 5 mJy, the lens- 
corrected normalized number counts are consistent, within error 



bars, with the GOODS-N counts (Berta et al. 



2010 1 showing 



evidence of a peak between 5-10 and 30mJy. Thanks to grav- 
itational lensing, we could extend the counts down to ~ 1 mJy 
at lOOjum and ~2mJy at 160yum, confirming the negative sub- 
Euclidean slope below the turnover. 

The Abell 2218 differential counts (not normalized to the 
Euclidean slope) at lOOyum show a faint-end slope similar to 
GOODS-N, but they are higher at fainter fluxes. At 160 /im the 
differential counts show a steeper slope than GOODS-N (- 1.82 
instead of -1.67 assuming the functional form of the counts: 
dN Ids cc S"). This is also reflected in the Euclidean-normalized 
counts. The high error bars seen at the faint end come from the 
poor statistics. The counts could be slightly underestimated at 
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Fig. 3. Number counts at 100 and 160 //m with lensing correction (red filled circles), normaUzed to the Euclidean slope, against the 
prediction of various synthetic counts models. See references in'Berta et al. ( |2010 l. Errors refer to pure Poisson statistics at 68% 
confidence limit GOODS-N counts from Berta et al. are contrasted in the shaded area. 



Table 2. Resolved CIB surface brightness for Abell 2218, and 
combined with the COSMOS contribution in the higher flux 
range ("PEP"). 





flux-range 
(mJy) 


(nW 


CIB 

mr^ sr"') 


d(CIB) 


fraction 

(%) 


100yumA2218 
lOOyumPEP 


0.94-35 
0.94-142 




7.9 
8.8 


3.4 
3.5 


55 ±24 
62 ±25 


160yumA2218 
leOyumPEP 


1.47-36 
1.47-179 




9.9 

11.3 


3.9 
4.1 


77 ±31 
88 ±32 



160yum at the lowest fluxes because the source density is ap- 
proaching confusion with no correction made for this. 



The counts are reproduced well by the models of Valiante 
[eraL] ( |2009l ) or |Rowan-Robinson et al.] ( |2009l ), both at 100 and 
160/zm, but other backwards evolutionary models similarly re- 
produce the downward turn below lOmJy. The summed contri- 
bution of resolved galaxies provides a lower limit to the IR back- 
ground and can be compared to the estimation of the CIB. Here 
we have adopted the l atest measure of it s surface brightness from 
COBE/DIRBE maps (" Dole et al.|2606| : 14.4 + 6.3 nWm^ sr 



at 100 jum and 12.75 ± 6.4 nWm~^sr~ at 160 //m (which is an 
interpolated value). 

The contribution to the CIB by Abell 2218 background 
galaxies above the 4cr detection threshold was computed by 
simply integrating the observed number counts (Table|2]i. Errors 
were computed by integrating the envelope of the counts with 
their uncertainties. More than half of the DIRBE CIB ( jDole et| 
[aL]|2006| has been directly resolved, which is also consistent 
with the surface brightness found by Bethermin et al. (20101 in 
a stacking analysis of 5/?/rzer/MIPS sources. 

Combining the deep counts in Abell 2218 with the results ob- 
tained by Berta et al (2010) in wider/shallower PEP fields (e.g. 
COSMOS), we were able to extend the integration flux range to 
142 mJy at 100 /vm and 179 mJy at 160 /im. Consequently, the 
resolved CIB fractions increase to 62 ± 25% and 88 + 32% in the 
two bands. One must keep in mind that not only is the CIB sur- 
face brightness from PEP affected by large uncertainties (cosmic 
variance), but the reference values by Dole et al. (2006 1 are also 
defined only within a factor of ~2. 

Another 9 massive lensing clusters will be targeted as part 
of PEP and another 40 more in the open-time Herschel Lensing 



Survey KP ( Egami et al.||20I0 l. We expect that, in the coming 
years, combined results from many lensing clusters will greatly 
improve the statistics of highly amplified sources and constrain 
source densities around ~1 mJy or below. By penetrating below 
the unlensed Herschel confusion limit and probing the high- 
redshift galaxy populations beyond the sensitivity limit of blank- 
field surveys, the fraction of the resolved CIB will increase. 
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